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models of interacting particles. In physics, the kinetic equations allow Lyapunov functionals (en¬ 
tropy, free energy, etc.). This may be considered as a sort of inheritance of the Lyapunov func¬ 
tionals from the microscopic master equations. We study nonlinear Markov processes that inherit 
thermodynamic properties from the microscopic linear Markov processes. We develop the ther¬ 
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1. Introduction 

1.1. What is the proper nonlinear generalization of the Markov processes? 

First order kinetics (the Kolmogorov-Chapman or master equation) is used in defining of nonlinear 
kinetic equations: the microscopic dynamics is replaced by Markov processes and then the large 
linear system is reduced to a nonlinear kinetics of some moments with referring to the law of large 
numbers for the stochastic evolution. This approach became very popular after the works of Kac 
ETll and Prigogine and Balescu ll42ll . In this sense, the Markov processes serve as a source of 
nonlinear kinetics. The stochastic simulation of chemical reactions ITT^ made the master equation 
approach to kinetics more popular in many applications. 

At the same time, master equation is considered as a simplest kinetic equation because it typi¬ 
cally defines a contraction semigroup. If we consider, for example, Markov transitions between a 
finite number of states, Ai —>■ Aj, then the probability distribution relaxes exponentially to an equi¬ 
librium and if the digraph of transitions is connected then the normalized equilibrium distribution 
is unique. On contrary, the interaction between states may produce nonlinear kinetic equations 
with various non-trivial dynamic effects. For example, if we write for two states (‘rabbits’ and 
‘foxes’) rabbit —)■ 2rabbits, fox -f rabbit —)■ (1 -f a) foxes (the interaction step), and fox —> 0, 
and apply the standard mass action law then we get the predator-pray Lotka-Volterra system with 
oscillations. 

The classical mass action law (MAL) systems are dense among the differential equations which 
preserve positivity (different versions of this theorem are proven in If34l [TTlI . see also discussion 
in [[351). Therefore, if we aim to consider a general class of kinetic equations which includes the 
MAL systems then the only important restriction is preservation of positivity. On this way we 
approach the theory of nonlinear Markov processes ifTlI^. 

In general spaces of states, the nonlinear Markov processes are the measure-valued dynamical 
systems which preserve positivity. They can be represented as the law of large numbers limits 
of general Markov models of interacting particles. The sensitivity analysis for these nonlinear 
evolution equations, that is the systematic study of the smooth dependence on the initial conditions 
and other parameters via the study of linearized system around a solution were performed in ll2^ 
I3D1331. 

Linear Markov chains have many Lyapunov functionals. For a finite chain with equilibrium 
distribution P* = (p*) they have the form 

Hk{p\\n = . (1.1) 

where P = (pi) is the current distribution and h is an arbitrary convex function on the positive 
semi-axis. These functionals were discovered by Renyi in 1960 lf43]l and studied further by Csiszar 
lO, Morimoto [[3^ and many other authors (see review in lU^ ). The functions Hh{P{t)\\P*) 
monotonically decrease (non-increasing) with time on the solutions P{t) of the corresponding 
master equations. Proposition [21 of Appendix extends the Morimoto result to continuous state 
models. 
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In physics, the kinetic equations allow Lyapunov functionals (entropy, free energy, etc.). This 
may be considered as a sort of inheritance of the Lyapunov functionals from the microscopic 
master equations. In this paper, we study nonlinear Markov processes that inherit thermodynamic 
properties from the microscopic linear Markov processes. We develop the thermodynamics of 
nonlinear Markov processes and analyze the asymptotic assumption, which are sufficient for this 
inheritance. 

1.2. Preliminaries: MAL, detailed balance and i7-theorems 

The classical thermodynamics follows the Clausius laws [|71 

1. The energy of the Universe is constant. 

2. The entropy of the Universe tends to a maximum. 

In practice, we assume that the ‘Universe’ is the minimal system, which is isolated with acceptable 
precision and includes the system of interest. 

Kinetics is expected be concordant with the laws of thermodynamics. In physical kinetics, 
Boltzmann’s H theorem established a link between the statistical entropy of one-particle distri¬ 
bution function in gas kinetics and the thermodynamic entropy [|3l. Boltzmann’s proof of his 
LT-theorem used the principle of detailed balance: At equilibrium, each collision is equilibrated by 
the reverse collision. This principle is based on the microscopic reversibility, the Newton equa¬ 
tion of motion for particles are invariant with respect to a time reversal and a the space inversion 
transformations. Five years before Boltzmann, Maxwell considered detailed balance as a conse¬ 
quence of the principle of sufficient reason ||T7]| . Later on, this principle was declared as a new 
fundamental law (Shi. For modern proofs and refutations of detailed balance we refer to |[T4ll . 

After Boltzmann, new kinetic equations in physics are always to be tested for concordance with 
the laws of thermodynamics. Many particular LT-theorems have been proved for various classes 
of kinetic equations. The principle of detailed balance has been widely used in these proofs. For 
MAL with detailed balance, the iT-function and the entropy production formula are very similar to 
the Boltzmann equation with detailed balance. Let Ai,..., be the components. For any set of 
non-negative numbers a pi, 13 pi > 0 (f = 1,... ,n, p = 1,..., m) a reversible reaction mechanism 
is given by the system of formal equations: 

OpiAi -f ... -f (VpnAn -i f3piAi (3 pn-^n • (1.2) 

According to the principle of detailed balance, each reaction has an inverse one and we join them 
in one reversible reaction. A non-negative real variable, concentration Ci, is associated with each 
component Aj, two positive constants, rate constants are associated with each elementary 
reaction and reaction rates are defined as 



2=1 2=1 
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The reaction kinetics MAL equations are 

| = (1.4) 

p 

where c is the vector of concentrations with coordinates c* and 7p is the stoichiometric vector of 
the elementary reaction, yp, = (3pi — api (gain minus loss). 

The principle of detailed balance for the MAL kinetics means that > 0 and there exists a 
positive point of detailed balance c* (c* > 0), where 


rp(c*) = rp(c*}{ = 


) (= r*), i.e. k+ = kp 


= r 


(1.5) 


2=1 


2=1 


For a given positive point of detailed balance, c*, the reaction rates include m independent positive 
constants, equilibrium fluxes r*, instead of 2m rate constants fcp : 


r+(c) = 


)=’-;n 


r- \ / n \ 

\ _ / \ ■4r T r / k^2 


, [C) = 




c 


( 1 . 6 ) 


i=l \ 2 - / 2=1 

/^-theorem for MAL kinetics with detailed balance is similar to Boltzmann’s LT-theorem. Take 


H{c) = Q ( In 


1 . 


(1.7) 


Simple calculation gives that for the kinetic equations (11.41) with reaction rate functions (11.61) 

-r;(c)) (In r+(c) - In r;(c)) < 0 (1.8) 


dH 

dt 


and 

dH 

= 0 if and only if (c) = r~ (c) for all p 

because (x — y) (In x — In y) > 0 for all positive x, y and it is zero if and only ifx = y. Hence, if 
there exists a positive point of detailed balance than H(c) decreases monotonically in time and all 
the equilibria are the points of detailed balance lf46]l . 

Physically, the constructed equations correspond to chemical reactions in a system with con¬ 
stant volume and temperature. For other classical conditions (isobaric systems, isolated systems, 
etc, the Lyapunov functionals are also known (see, for example, [|26ll48]| ). 

For many real systems the reaction mechanism includes both reversible and irreversible reac¬ 
tions. For them some reverse reactions are absent in the reaction mechanism (11.21) . (It is convenient 
to use such notations that all direct reactions are present and some reverse reactions are absent). 
The systems with irreversible reactions which are the limits of the fully reversible systems with 
detailed balance when some of the equilibrium concentrations tend to zero are described [l20ll^ . 
If the reversible systems obey the principle of detailed balance then the limit system with some 
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irreversible reaetions must satisfy the extended principle of detailed balance. It is proven in the 
form of two conditions: (i) the reversible part satisfies the principle of detailed balance and (ii) the 
convex hull of the stoichiometric vectors of the irreversible reactions does not intersect the linear 
span of the stoichiometric vectors of the reversible reactions. These conditions imply the existence 
of the global Lyapunov functionals and alow an algebraic description of the limit behavior. The 
extended principle of detailed balance is closely related to the Grigoriev - Milman - Nash theory 
of binomial varieties [[24ll . 

1.3. Thermodynamics beyond detailed balance 

In the original form of iJ-theorem the microscopic reversibility (invariance of the microscopic 
description with respect to time reversal) is used to prove the macroscopic irreversibility, the exis¬ 
tence of the time arrow {H decreases monotonically due to kinetic equations). Elegant paradoxical 
form of this reasoning leaves, nevertheless, concern about its generality: does the macroscopic 
irreversibility need the microscopic reversibility? In 1887 Lorentz formulated this concern explic¬ 
itly. He stated that the collisions of polyatomic molecules are irreversible and, therefore, Boltz¬ 
mann’s LT-theorem is not applicable to the polyatomic media iHTI . Boltzmann found the solution 
immediately and invented what we call now semidetailed balance or cyclic balance or complex 
balance BUl. For the Boltzmann equation this new condition allows a nice schematic representation 
(see Figure [U for detailed balance and Figure [2] for complex balance). Now, it is proven that the 
Lorentz objections were wrong and the detailed balance conditions hold for polyatomic molecules 
m- Nevertheless, this discussion was seminal and stimulated Boltzmann to discover new general 
conditions of thermodynamic behavior. 



Figure 1: Schematic representation of detailed balance for collisions. The four-tail scheme repre¬ 
sents intensity of the equilibrium flux of collisions with given velocities. 



z 

inputs 



V 


W 


Figure 2: Boltzmann’s cyclic balance is a summarised detailed balance condition: at equilibrium 
the sum of intensities of collisions with a given input v + w —)■ ... coincides with the sum (or 
integral) of intensities of collisions with the same output... v + w. 

For the MAL kinetics, the Boltzmann cyclic balance condition was rediscovered in 1972 [|25]|. 
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It got the name “eomplex balanee” (balanee of eomplexes). The eomplex balanee eondition has 
the form of summarised detailed balanee (Figure O. 

Consider a reaetion meehanism in the form 


OipiAi + ... + apnAn —>■ f3piAi + pnAfi 


(1.9) 


Here, the reverse reaetions if they exist partieipate separately from the direet reaetions and re¬ 
versibility is not eompulsory. This form is eonvenient for systems without detailed balanee. The 
MAL reaetion rate is Vp = kp Hj and the kinetie equations have again the form (11.41) . 

A positive eoneentration veetor c* is an equilibrium if = 0’ where Vp = kp Hi(^1 

Complexes are the formal sums in the left and right hand sides of (11.91) . There are 2m veetors 
of eoeffieients ap = {upi) and j3p = {j3pi) (p = 1,... ,m). Some of them might eoineide. Let 
{yi ,..., yq} be the distinet eoeffieient veetors: for eaeh yj there exists sueh p that yj = ap or 
yj = l3p, and for eaeh p, ap there exists sueh j, I that yj = ap and yi = l3p. 

A positive point c* is a point of eomplex balanee if for eaeh yj 


Z = (1.10) 

p^yj=o‘p p,yj=/3p 

This is exaetly the summarized detailed balanee eondition (eompare it to Figure [2l). The eomplex 
balanee eonditions (11.101) are suffieient for the LT-theorem: dH/dt < 0. To demonstrate this 
inequality, we eonsider the deformed stoiehiometrie meehanism with the stoiehiometrie veetors 
whieh depend on parameter A G [0,1]: 

+ (1 ~ •^)/^p ) ^p{^) = -^/^p + (1 ~ A)Q;p. 


Introduee an auxiliary funetion 6{\), that is the sum of the reaetion rates of the deformed meeha¬ 
nism (with the same equilibrium fluxes). For a given eoneentration veetor c 


p i=l 



Simple ealeulation gives 


dH 

dt 


dX 


A=1 


The funetion 6*(A) is eonvex and the eomplex balanee eonditions (11.101) imply 6*(0) = 6^(1), there¬ 
fore under these eonditions ^'(1) > 0 and H monotonieally deereases in time. 

For first order kineties (eontinuous time Markov ehains or master equation) the eomplex bal¬ 
anee eonditions are just the stationarity eonditions (the so-ealled balanee equations) and hold at 
every positive equilibrium. This gives immediately the LT-theorem for first order kineties with 
positive equilibrium c* and without any additional eonditions. 
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1.4. From Markov kinetics to MAL with complex balance condition 

The semidetailed balance conditions (Figure[2l) for Boltzmann’s equation were produced by Stueck- 
elberg ||45]| from the Markov model of collisions (Stueckelberg used the S'-matrix notations and 
presented the balance equation as the unitarity condition). 

MAL for catalytic reactions with a priori unknown kinetic law was obtained in the famous work 
of Michaelis and Menten lIMll . They postulated that substrates form complexes (‘compounds’) 
with enzymes, these compounds are in equilibrium with enzymes (fast equilibria), and the con¬ 
centrations of the compounds is small. The compound-substrates equilibria can be described by 
equilibrium thermodynamics and the kinetics of the compounds transformations is just a Markov 
chain (because for very small concentrations of reagents only first order reactions survive). These 
asymptotic assumptions lead to MAL. Michaelis and Menten studied very simple reaction, there¬ 
fore the additional relations between reaction rate constants did not appear but the Stueckelberg 
approach extended to the general reaction kinetics gives the semidetailed balance (complex bal¬ 
ance) condition Il22ll . 

The asymptotic assumptions of the Michaelis-Menten-Stueckelberg works are illustrated by 
Figure[3l Each complex (^. apiAi or l3piAi) is associated with its compound (Bp in Figure). It 
is assumed that the complex is in fast equilibrium with its compound and the concentration of com¬ 
pounds can be found by conditional minimization of a thermodynamic potential (under isothermal 
isochoric condition this is the free energy). The second asymptotic assumption means that the con¬ 
centrations of compounds are small with respect to the concentration of reagents. This condition 
allows us to find the concentrations of different compounds independently. For the perfect free 
energy F = constant x H with H given by (11.71) these concentrations might be found explicitly. 

It should be stressed that the fast equilibrium assumption was later eliminated from enzyme 
kinetics by Briggs and Haldane O and what is often called the Michaelis-Menten kinetics is the 
Briggs-Haldane kinetics (for the modern analysis of this system we refer to Il44l f. Nevertheless, 
the idea of intermediate complexes which are in fast equilibria with stable reagents is crucially 
important for production of dynamic MAL from thermodynamics. This idea was reanimated and 
systematically used in the theory of the activated complex and reaction rates flU [TOl |40l . There¬ 
fore, the Michaelis-Menten-Stueckelberg limit (Figure [3]) may be called the Michaelis-Menten- 
Stueckelberg-Eyring limit. 

In our work, we study the Michaelis-Menten-Stueckelberg limit of Markov processes with 
general space of states and obtain for them the generalized MAL (GMAL) with thermodynamic 
properties. 


2. Derivation of MAL and GMAL on the arbitrary state space 

2.1. Thermodynamics of particles 

To set a scenery suppose a species or a particle can be represented by a point x in a locally compact 
metric space X with some fixed Radon measure M(dx). The distribution of (possibly infinitely 
many) particles in X can be specified by a finite measure. We shall deal only with distributions 
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Figure 3: The Miehaelis-Menten-Stueekelberg limit. 


that have densities (eoneentrations) c{x) G (M) with respeet to M. 

Let the thermodynamie properties of a eoneentration c be eharaeterized by the ’free energy’, 
whieh is given generally by a funetional F(c(.)) defined on L^{M) (or some its subspaee, the 
domain of F). We shall assume that F is smooth in the sense that the variational derivative 5F/5c 
(for positive c) exists with respeet to M, defined by the equation 


A 

dh 


F{c + huj) = 


/i=0 


IX 


6F 

6c{x] 


■u{x)M{dx) 


for c, CO from the domain of F. 

For many applieations, 5F/5c has a logarithmie singularity as c —?■ 0. Further on we assume 
the positivity of c when it is neeessary. 

Two basie examples that eover all known physieal models should be kept in mind. In the first 
one 

F(c(.)) = [ 'tl}{c{x))M{dx) (2.1) 

Jx 

with a funetion -0, whieh is smooth on positive arguments. In this ease 

=f(c{x)). 


6c{x] 


This ineludes the ease of finite X = {1, • • • ,k} with 


F(ci,--- ,Cfc) = 

j 


As another partieular ease let us mention the standard perfeet gas free energy given by 


m.)) 





6F 

6c{x) 



( 2 . 2 ) 


with some equilibrium distribution c* (we omit here the eonstant faetors). 
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Of interest is also its more general version, so-ealled ’generalized entropy’ function (with in¬ 
verted sign) 



(2.3) 


where h is any convex smooth function on R+ = {x > 0}. Function (12.21) is obtained from (12.31) 


for h{x) = x(lnx — 1). Choosing h{x) = — Inx leads to the so-called Burg relative entropy 



(2.4) 


In the second example X = x {1, • • • , /c} with M{dx) being Lebesgue measure on each 
component (more generally, instead of R*^ one can use a manifold, but we shall stick to R*^ for sim¬ 
plicity). With some abuse of notation we shall denote the elements of X by a pair {x, j), x G R'^, 
j = 1, • • • , /c or sometimes by Xj. The concentration c becomes a vector c = (ci(a;), • • • , Ck{x)) 
with its gradient Vc(a;) = (Vci(a;), • • ■ , Vcfc(a;)), where 



The free energy is specified by the equation 



(2.5) 


with some smooth function This includes the case of finite X = {1, • • • , fc} with 


F(c(.)) = 'i/;(ci,-- - ,Cfc). 


( 2 . 6 ) 


The well established particular case of (12.51) is 



(2.7) 


used in the Chan-Hilliard model of diffusion. The simplest version of (12.51) is the decomposable 
case: 



k 


( 2 . 8 ) 


i=i 

The variational derivative for F of type (12.51) is the standard Euler-Lagrange one: 

5F d dif) 



(2.9) 


Scj{x) dcj ^dx^dViCj 
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2.2. Compounds 

Our objective is to describe the process of transformation of particles (chemical reactions, colli¬ 
sions, etc). The collections of k particles can be given by points in SX^ and the collection of an 
arbitrary number of particles in SX = where SX’^ and SX are the quotient-spaces of 

and X = with respect to all permutations. Symmetrical probability laws on X^ (which are 

uniquely defined by their projections to SX^) are called exchangeable systems of k particles. With 
some abuse of notations we shall use the same bold face letter notations, say x = (xi, • • • , Xk), 
both to denote the points of X^ and S'X^. We shall also use the notation X>k = Uj>fcXT 

The main idea of the intermediate state or the activation complex assumption (that we shall 
adopt here) is that any reaction changing the collection of k particles x = (xi, • • • ,Xfc) to the 
collection of / particles y = (|/i, • • • , |/fc) is not a one step operation, but the three step one: before 
the interaction is enabled the k particle x = (xi, • • • , x^) should form the intermediate state x, 
which we shall call a compound of size k (consisting of the same k particles x), then the compound 
X turns to the compound y, which in turn can be dissolved into its components y: 


X = (xi • • • , Xfc) X y y = (l/i, • • • , l/i). 


( 2 . 10 ) 


This concept of the intermediate states allows one to speak about the distribution of compounds 
present in the system. Introducing some fixed symmetric measures on SX^ allows one to 
reduce attention to distributions specified by the densities (concentrations), which, for compounds 
of size k, are given by the symmetric functions Cfc(x) G L^{Mk). 

Let us denote by M the measure on X with the coordinates (M = Mi, M 2 , • • •) and by 
p{x,dy) the stochastic kernel on X with the coordinates pk{x,dy) so that for a function / = 


(/i,/2,---) onX, 



The supports of measures Mj specify the set of compounds that can be formed in the system. 

Of course, the simplest example of the measures Mk are the projections on SX^ of the products 
M{dx) 0 • • • 0 M{dx) (k times), but this example is not sufficient (as we shall see below) to cover 
all cases of interest. However, in order to develop a theory, some link between M^ and M should 
be made. Our main assumption about M^ will be that the projection of all Mk on the one-particle 
states is absolutely continuous with respect to M, namely 



( 2 . 11 ) 


with a symmetric stochastic kernel /ifc(x, dx), x G X, x G SX^ By symmetry, (12.1 II) rewrites 
as 



( 2 . 12 ) 


for any j. 
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This assumption is crucial for the possibility to relate the concentration of partieles with the 
eoneentration of eompounds. Namely, as shown in Appendix (see Proposition [3]), the kernels /i 
in (12.1 II) can be ehosen in sueh a way that if • • • , Xk) = Cfc(x) is the eoneentration of the 

eompounds x of size k, the eoneentration of partieles involved in these eompounds equals 

/ Cfc(a;, X2,--- , Xk)fik{x, dx2--- dxk). (2.13) 

Notice that the coefficient 1/k was introdueed in (12.1 II) to avoid any coefficients in (12.131) . The 
kernels ^k appearing in Proposition |3] will be ealled stoichiometric kernels, as they present natural 
analogs of the stoiehiometrie eoeffieients of the theory of ehemieal reaetions on a finite state spaee. 

Therefore, if ({x) is the eoneentration of free partieles (not involved in the eompounds), the 
total eoneentration of partieles is 

OO n n 

c{x) = C{x) + '^ C{x,X 2 ,-■ ■ ,Xk)tik{x,dx 2 -■-dxk) = c{x) + / C{x,y)ti{x,dy), 

k=2 -JSX 

(2.14) 

where C(x) = Cfc(x) for x E and /i(x, dy) is the stoehastie kernel on X (stoiehiometrie kernel) 
with the ’eoordinates’ t^k{x, dy) on SX^. 

Further on we shall assume for simplieity that the size of possible eompounds is uniformly 
bounded, so that all sums over sizes used below are finite. This restrietion is also natural from the 
praetieal point of view, as the sizes of eompounds met in practice are very small (usually 2 or at 
most, and rarely, 3). 


2.3. QE and QSS 

The quasi-steady-state (QSS) assumption states that the compounds exist in very small concen¬ 
trations as eompared with the concentration of free particle (beeause they form and dissolve very 
quiekly) and the quasi-equilibrium (QE) assumption states that the reaetion of equilibration be¬ 
tween partieles and eompounds is mueh faster than the reaetion between eompounds meaning that 
the eompounds exist all the time in a fast equilibrium with the set of basie partieles. Let us diseuss 
the important eonelusions from these assumptions. 

First of all, by QSS, the free energy of the eompounds ean be taken in the form of the perfeet 
free energy (the free energy of the ideal gas or of dilute solutions), so that the total free energy of 
the system beeomes 

Ftot{c, C 2 , C3, • • •) = = nc{-)) + f C(x) fin ^ - l) M(dx) 

JSX:>2 V S (XJ / 




Mj{dx) 


(2.15) 


with F the free energy of partieles as introdueed above and with Q, j > 1, some equilibrium 
eoneentrations. Generally speaking, should depend on c, but again by QSS, the concentration of 
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particles are large and vary slowly as compared with the eompounds implying that this dependence 
can be neglected in the first approximation. In the same approximation we do not distinguish the 
concentrations of free particles Ci(^) their total concentration c(x). 

By QE the eompounds are all the time in equilibrium with partieles. As in equilibrium the free 
energy takes its minimum, Q, j > 2, ean be found from the eondition of the extremum: 


^ (c(-)(C + ea;)(y)) (^In M(dy) : 

that should hold for all symmetrie functions a;(x) on A’> 2 . By the definition of the variational 
derivative this implies 



0 , 


'X 


6F 

6c{x) 


I sx 


uj{x, y)fi{x, dy)M{dx) + 


f ‘^Wto|^M(dx) = 0 

IsXy, C*(x) 


(2.16) 


for all u and hence, by (12.111) . 


U). 


{xi,X2, ■■■ , Xj) l^ln Mj{dxidx2 ■ ■ ■ dxj) = 0 


(2.17) 


>SXi V 

for all j > 2 and Uj. By the symmetry this implies 




(2.18) 


Consequently 


In 


Ci(x) 




6F 


c;(x) ^6c{xi)' 

so that finally, for any j > 1, the minimizing concentrations are 


(2.19) 


0(x;c)=C;(x)exp||:^|, 


( 2 . 20 ) 


2.4. Dynamics 

The next eonsequenee of the QSS is that the dynamics of compounds should be linear, because, 
their eoneentration being small, one ean negleet their interaction. This includes the dynamies of 
free partieles, as their interaction has been aeeounted for by the formation of eompounds. 

By (17.101) . assuming (17.31) and using the notations for the concentrations of compounds intro¬ 
duced above, the general jump-type Markov evolution on eoneentrations ( ean be written in the 
eoneise form 
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C(x) = / [C(y)z^(x,rfy)-C(x)i/(x,rfy)]. (2.21) 

JSX 

Here u is some eolleetion of transition kernels 

z/(x,c(y) = {i/fc^z(x,c(y), x G SX^, y G SX^} 

and the eolleetion i) is defined via the following equality of measures on eaeh pair SX^ x SX’': 

z/fc^z(y, dx)Mfc(dy) = dy)Mj(dx). (2.22) 

The additional eonstraint arising from thermodynamies eomes from our assumption, see (12.151) . 
that (* are equilibrium eoneentrations of eompounds and therefore they should supply equilibrium 
for their linear evolution (12.211) . that is 


/ [C*(y)z>(x, dy) - C*(x)z/(x, dy)] = 0 (2.23) 

JSX 

for any x. 

It is useful to distinguish a subclass of processes where particles themselves do not serve as 
compounds, or in other words, direct transitions X —)■ SX >2 and SX >2 —)■ X are not allowed: 


i/(x, dy) = 0, i/(y, dx) = 0, x e X,y e SX> 2 - (2.24) 

If this is the case condition (12.231) should be understood as 

[ [C*(y)^(x, dy) - C(x)z/(x, dy)] = 0, x G SX> 2 . (2.25) 

JSX>2 

Otherwise, for (12.231) to make sense, equilibrium quantities C*{^)^ x E X, should be defined 
somehow to complement the definitions of C*(x; c) for x G SX >2 given by (12.191) . 

A simpler subclass of processes worth being mentioned present the evolutions preserving the 
number of particles in the compounds. In this case, the evolution (12.211) decomposes into the 
independent evolutions in each SX^: 


4(x) = / [Cfc(y)z>fc^fc(x,dy) - Cfc(x)i^fc^fc(x,dy)], /c = l,2, •••. (2.26) 

Jsx*^ 

By (12.141) . the evolution of the compounds (12.211) implies the following evolution of the total 
concentration c: 


c{x) = / [C(y)z>(a:, dy) - C{x)y{x, dy)] 


’sx 


+ / fx{x,dy) [C(z)j>(x,y,dz) - C(x,y)z/(a:,y,dz)]. 


’SX 


'SX 


(2.27) 
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It remains now to put it all together. As shown above, by the QE assumptions Cfc for k > 1 are 
expressed by (12.191) in terms of c. By QSS, for fc = 1 we have approximately 

Cl ( 3 ^) = Cl ( 3 ^; c) = c(x). (2.28) 

Finally, apart from the transformation of partieles it is natural to allow additionally their movement 
in X aeeording to some Markov proeess with the generator L (only free partieles are moving, as 
the movement of the short-lived eompounds ean be negleeted). Then the final evolution of the 
eoneentration beeomes 


c{x) =L*c{x) + / [C(y; c)i){x, dy) - Ci(x; c)u{x, dy)] 


I sx 


(2.29) 


'sx 


f^ix,dy) / [C(z;c)z>(x,y,dz) - C(x,y;c)z/(x,y,dz)] 


isx 


supplemented by (12.191) and (12.281) . This evolution ean be ealled the generalized mass action law 
(GMAL). It is the extension to an arbitrary state spaee X of the finite-state-spaee GMAL. The latter 
was developed in general by Gorban et al [fT3l [TTl |22l following the ideas of Miehaelis-Menten, 
Eyring, Stueekelberg, and many others. For the diffusion equations the formalism of GMAE was 
also elaborated [|2T]| . 

If eondition (12.241) holds (partieles are not eompounds), evolution (12.291) rewrites in a simpler 
form 


c{x) = L*c{x)+ fi{x,dy) [({z; c)i>{x,y, dz) - C{x,y, c)iy{x,y, dz)]. (2.30) 


'sx 


'sx^ 


If the partieles themselves are given in small eoneentrations, so that their free energy has the 
perfeet form (12.21) . evolution (12.291) turns to 


c(x) =L*c(x) 


J sx . 
ti{x,dy) 


'sx 


dy) - c(a;)^^z/(a:, dy) 
c*(y) c*(a;) 

c(z)^;^z>(a;, y, dz) - c{x, y) ^}^'^l u{x, y, dz) 
'sx L ^ ^ y) 


(2.31) 


where C*{x) = c*{x) for x G X and 


c(xi, ■■■ ,Xk) = c(xi) • • • c(Xfc), c*(xi, • • • , Xk) = c*{xi) ■ ■ ■ C*{Xk). 

This is the evolution of the MAE for an arbitrary state spaee X. 

Important to observe that, for the MAE evolution (12.311) . equilibrium quantities C*{x) = c*{x) 
for X G X are explieitly speeified from the expression of free energy, and thus the eondition 
(12.231) is well defined without the restrietion (12.241) . Moreover, eondition (12.231) supplemented by 
the similar eondition on the free evolution of c, that is assuming L*c* = 0, implies that c* are 
equilibrium eoneentrations to (12.311) . 
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2.5. Basic examples 

For the case of a discrete state space X = {Ai , • • • , it is more convenient (and well established 
in the literature) to use separate enumeration of compounds and reactions. For each reaction r, 
denoting by the number of particles A* entering the input compound B~ and by the number 
of particles Aj entering the output compound the reaction can be described schematically as 

'y ^ OlriAi ^ —)■ ^ y ^ (3riAi. 


Here ari,Pri are known as stoichiometric coefficients and the vector Ur = {Ph — c^ri) as the 
stoichiometric vector of the reaction r. In this notation the evolution (12.291) becomes 


f-i y y ^jiCi (2.32) 

¥j 

with some Kij playing the role of transitions /i(x, dy) of (12.291) : the summation is over all pairs of 
compounds (/, j) and 

(c)z/z,| (2.33) 

for each compound /, with F, xjj from (12.61) . 

Most of real life evolutions involve compounds consisting of only two particles. If only pairs 
to pairs transitions can occur then GMAL (12.291) and MAL (12.311) take the form 



c{x) = L*c{x)+ I ii{x,dy) / [({zi, Z 2 ]c)v{x,y,dzidz 2 ) - C{x,y]c)v{x,y,dzidz 2 )], (2.34) 

Jsx^ 


' X 


and respectively 

c(x) = L*c(x) 


y{x,dy) 


lx 


c(zi)c(z2) , , N c(a;)c(?/) ,, , , , , ^ 

Z2)y[x, y, dzidz2) - , . , . ({x, y)u{x, y, dzidz2) 


sx^ [c*{zi)c*{z2) 


c*{x)c*{y) 


(2.35) 


Notice now that the densities entering the kernel y can be transferred to the rates v. Hence, as 
was already pointed out, basically only the support of y is essential. It turns out that two particular 
cases cover all interesting examples. The first comes from the assumption that any pair of particles 
can interact. In this case, the measure M 2 on pairs can be taken to be proportional to the product 
measure M{dx)M{dy) and then one can take 

y{x,dy) = M{dy). (2.36) 

In the second case, the state space X is the product R"* x V equipped with the measure dxM{dv), 
where the first component is interpreted as the position in space and where it is assumed that a pair 
of particles can interact if and only if their positions in space coincide. In this case one has 

y{{y, v), d{z, w)) = S{y - z)M{dw). (2.37) 
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For instance, the full Boltzmann equation is of that type. 

Let us also distinguish two oases of interest oonoerning transitions u. The simplest ease is of 
oourse when the transitions i'{xi,X 2 ', dyidy 2 ) are absolutely oontinuous with respeot to the produot 
measure M{dyi)M{dy 2 ). If also (12.361) holds evolution (12 341) turns to 

c{x) = L*c{x)+ M{dy) / [({zi, Z 2 ; c)u{zi, Z 2 ; x,y)-Cix,y; c)i^{x,y; zi, Z 2 )]M{dzi)M(dz 2 ) 

Jx Jsx'^ 

(2.38) 

This example is however of rather limited applioability. More interesting situation ooours when 
there is given another measure spaoe with the measure du and a family of M 2 -measure-preserving 
bijeotions : X'^ —)■ depending on a; as a parameter suoh that 

u{x,dy)= f B{x,uj)5{y - Gui{x))dudy (2.39) 

Jn 

with some funotion B{x, ui) on x 12. Sinoe 


/ / f{y,x) 5 {x-G^{y))B{y,u)dxduM 2 {dy) = 

Ix^ Jn 

f{GZ\^)^^)B{G~\x),u)duM2{dx.) = 

'x Jn 

it follows (see (17.31) ) that 



X Jn 



'x2 Jn 


/(y, G^{y))B{y, uj)duM 2 {dy) 
f{y,^)d{y-G~\x))B{y,u)dyduM 2 {dx.), 


z/(x,dy)= [ B{y,uj)6{y-Gj{x))dujdy. 

Jn 

Consequently, assuming again (12.361) . evolution (12.341) turns to 

c{x) = L*c{x)+ j M{dx 2 ) Ac((3^J^(x);c)5((:7j^(x),a;) - C(x;c)5(x,a;)]cicu. (2.40) 

Jx Jn 

In partieular, if B{x, oj) is invariant under the aetion of G^}, this simplifies to 

c{x) = L*c{x) + f M{dx 2 ) f B{x,uj)[({G~^{x.);c) - ({x;c)]duj. (2.41) 

Jx Jn 

For instanee, the spatially homogeneous Boltzmann equation is of that type, as well as the mollified 
Boltzmann equation and their fcth order extension, see |[29l . 

To give an example of evolutions arising from (12.371) . assume the rates u are absolutely eontin- 
uous with respeet to the seeond variable and are independent of the first variable. Then evolution 
(12.341) turns to 

c{yi,vi) =L*c{yi,Vi) + [ [ M{dv2)5{y2- yi) [ M{dwi)M{dw 2 ) 

7Rd Jv Jsv^ (2.42) 

X [C(w^i, yuW 2 , f/ 2 ; c)u{wi,W 2 ; Vi,V 2 ) - C(wi, yi] V 2 , 2 / 2 ; c)z/(ui, V 2 ; Wi, W 2 )]. 
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3. Analysis of equilibria 


3.1. Evolution of the free energy 

We are interested in eonditions ensuring the deerease of the free energy F(c(.)). If c(x) evolves 
aeeording to (12.301) . the free energy evolves as (where (12.111) is used to get rid of /r) 


/ kM,(dy)-^ / 

^ Jsxk dc{yi) Jsx >2 


[C(z; c)j>(y; dz) - C(y; c)u{y, dz)]M{dz). 


Using symmetry and introdueing a handy speeial notation 


( SF Y 


6F 


Uc(y)/ 

this rewrites as 

^(c(')) = f 4^L-c{v)M{dy) 


, y = ( 2 / 1 ,• • • , 1 /i) e x', ( = 1,2,•••, 




( [(;(z; c)v(y, dz,) - C(y; c)v(y, dz)] 


or, using the definition of u, as 

F(c(.))^j^,^L>c(y)M(dy) 


+ 


(dc^) dz)M(dy)]. 


JSX>2 JSX>2 

Finally, relabeling the variables in the seeond term of the last integral yields 

6F 


F{c{-)) = / ^L*cix)Midx) 

Jx oc{x) 


6F 


I SXy2 JSX>2 I \Sc{y)J \6c{z 

Turning to the general ease (12.291) we find similarly that 

5F 


(3.1) 


C(z;c)i/(z,dy)M(dz). 


F{c{.)) = / ^L*c{x)M{dx) 
Jx oc{x) 


isx Jsx 


(—)' 

vMy)/ 


6F 

6c{z] 


(3.2) 


C(z;c)i/(z,dy)M(dz). 
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3.2. Complex balance and detailed balance 

Evolutions (13.11) or (13.21) can be considered as continuous-state-space analogs of the discrete state- 
space representation giving the dynamics of the free energy in terms of the sum over reactions, 
as here we have the representation in terms of the integral over the pairs (y, z) that effectively 
parametrized possible reactions. 

With this analogy in mind, and dealing again first with evolution (12.301) and (13.11) . we can now 
generalize the trick used for the discrete case and introduce the auxiliary function 


6)(A) = 6»(A; c) = / / M(ciz)C*(z)z/(z; dy) 

JSX^2 Jsx>2 


(3.3) 


so that 


«'(A) = ^«(A;c) = 


'SX>2 Jsx>2 


M{dz)C{z)u{z;dy) 


6F 

6c(z) 




»"(A) = ^«(A;c) = 


’SX>2 JSX>2 


5c{z) 

M(dz)C*(z)iA(z;dy) 


X exp < A 


6F 

6c{z] 


+ (1 - A) 


6F 

6c{z) 

6F Y 


<^c(y)/ 


vMy)/ 


(—V 

vMy)/ 


(3.4) 


Hence 6*"(A) > 0, so that 9{\) is a convex function, and moreover, by (13.131) . 

F{c{.))=f ■^^L*c{x)M{dx) 

Jx oc{x) 


(3.5) 


Consequently, the conditions 



6F 

6c{x) 


L*c{x)M{dx) < 0 


(3.6) 


and 

d(0) < d(l) 


(3.7) 


are sufficient for the decrease of free energy along the evolution (12.381) : F(c(.)) < 0. 

Inequality (13.71) introduced in ifTSll is referred to as (3^-inequality. It is a natural weakening of a 
stronger condition 

d(0) = d(l), (3.8) 
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which is often easier to analyze, since it rewrites as 

f f M(rfz)C*(z)z/(z;rfy)expi 

Jsx^2Jsx^2 [Voc(z^ 


or equivalently, again using the definition of i>, as 

J M(dz)[C*(z)z/(z;dy) - C*(y)f>(z;dy)]exp I 


= 0. (3.9) 


Assuming the functional F is rich enough, so that the linear eombinations of the exponents 

for all eontinuous funetions c are dense in the spaee of eontinuous funetions on SX> 2 , as is the 
ease for the MAL evolution, (13.91) implies 


'SX>2 


[C*(z)i^(z; dy) - C*(y)f>(z; dy)] = 0 


(3.10) 


for all z G SA’> 2 , whieh is the equilibrium eondition (12.251) . 

This eondition (13.101) is ealled the complex balance condition for evolution (12.301) . As was 
shown, together with (13.61) . it is suffieient for evolution (12.301) to ’respeet’ thermodynamies: F{c{.)) < 
0 . 

In partieular, for the pair-interaction dynamics (12.381) and (12.401) . the complex balanee condition 
takes the forms 




[C*(y)z/(y; x) - C*(x)i 2 (x; y)] dy = 0, x G 


and respeetively 


[C{y)B{y, w) - C*(x)5(x; u;)] dw = 0, x G 


(3.11) 


(3.12) 


and the evolution of the free energy 

F{c{-)) = [ ■^^L*c{x)M{dx) 
Jx oc{x) 


f f 

■ 5F 

SF 

1 

SF 

SF ' 

/x2 Jx '2 

.^c{yi) 

Sc{y2) 

Sc{xi) 

Sc{x2)_ 


C(x; c)i/(x; y)M(dx)M(dy) 

(3.13) 
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and respectively 




' X 


6F 

6c{x] 


L*c{x)M{dx) 


1 

4 


f [ 

■ 6F 

5F 

1 

5F 

5F ' 

/x2 7x2 

.^c{yi) 

< 5 c(|/ 2 ) 

6c{xi) 

5 c{x 2 ) _ 


C(x; c)5(x, uj)M{dx.)du, 

(3.14) 


where 

y = {yi,y2) = 

Of course (13.101) holds if 


C(y)Ky;x)-C(x)Kx;y) = o 


(3.15) 


for all X, y G SX> 2 - This more restrictive condition is called the detailed balance condition for 

(12301) . 

Turning to more general evolution (12.291) . (13.21) we shall reduce our attention only to MAL 
evolution, where 

cW = c-Wexp|^ 

for the free energy in the perfect form (12.21) . This makes the notations for ^(x) consistent with the 
notations Ci(2^) = c{x). Consequently, introducing 9 by the equation 


9{X) = 


M{dz)(*{z)u{z; dy) exp < A 


'sx J sx 


6F 

6c{z] 


+ (1 - A) 


6F Y 

My)/ 


M(dz)C(z)z/(z;ciy) 


'SX Jsx 


’c(z) ■ 

A 

1- 

1- 

o 

* 


_c*(y)_ 


l-A 


(3.16) 


yields again (13.51) . Moreover, condition 0(0) = 0(1) becomes equivalent to condition (12.231) . which 
is the complex balance condition for general MAL. 


3.3. Points of equilibrium 


Let us start with the MAL dynamics. As we know already, then c* is an equilibrium point. Are 
there other (positive) equilibrium points? Assume the complex balance condition (12.231) holds, 
and let c(x) be an equilibrium. Then we have 0(0) = 0(1) and 0'(1) = 0, which together with 
convexity of 0 implies that 0(A) is a constant (for given c). Hence 0"(A) = 0, and consequently, by 
(EH), 


SF 

Sc(z} 


(TL) =0 

vMy)/ 


(3.17) 


on the support of the measure M(dz)(^*(z)z/(z; dy), which coincides with the support of the mea¬ 
sure M{dz)i'{z] dy) if all (* are strictly positive. 
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In particular, it implies the following. Suppose the complex balanee condition (12.231) holds 
for a MAL evolution, all (* are strietly positive, the evolution preserves the number of particles 
and the measure u{z, dy)Mk{dz) on x has the full support for at least one k > 1. Then 
c is an equilibrium if and only if 5F/5c{x) is a eonstant (as a funetion of x) and L*c = 0, and 
henee, buy the strueture of F, if and only if c(x) eoineides with c* up to a multiplieative eonstant. 
Alternatively, assume (* > 0, (12.231) holds and the measure z/(z, dy)Mk{dz) on X^ x X^ has the 
full support for at least one pair k ^ 1. Then c* is the only (positive) equilibrium. 

Turning to evolution (12.301) and (13.11) we ean eonelude similarly that if eomplex balanee eondi- 
tion (13.101) holds, all C*(z) for z G SX are strietly positive, the evolution preserves the number of 
partieles and the measure i/(z, dy)Mk{dz), z, ?/ G has the full support for at least one k > 1, 
then c is an equilibrium if and only if 5F/5c{x) is a constant (as a function of x) and L*c = 0. 
Alternatively, assume Q* > 0, (13.101) holds and the measure z/(z, dy)Mk{dz) on X^ x has the 
full support for at least one pair k ^ 1. Then c is equilibrium if and only if 5F/5c{x) = 0 and 
L*c = 0. 

3.4. Comments on the transformations of the free energy 

The GMAL evolution (12.291) will not be ehanged if we make a linear shift of the free energy 
changing F to 

F = F + j 

with some uj{x) and simultaneously ehange C* to 

C(x) = C(x)exp | - , X = (xi, ■■■ ,Xj) e XT 

Redueing our attention for simplieity to evolution (12.301) . (12.241) . suppose the complex balance 
condition (12.251) does not hold. The natural question arises whether we ean find a funetion oj sueh 
that for new F, C,* it becomes valid, that is 



(y)}j>(x, dy) - C*(x) exp{-a;®(y)}z/(x, dy)] = 0, 


xG^X>2, (3.18) 


where 

a;®(x) = uj{xi) H-h uj{xk), x = {xi, ■■■ ,Xk). 

In diserete setting this question ean be effeetively answered algebraieally by the so-ealled de- 
fieiency zero theorem IfTTl In our setting let us note only that, if all initial C,* had a product 
form, this question reduces to the directly verifiable question on whether the produets 


C*(x) = c*(a;i) • ■■c*{xk) 


satisfy (12.251) . 
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4. More general free energy for compounds 


Motivated by Morimoto’s Theorem [2l ll^ . it is natural to use for the free energy of the compounds 
in (12.151) the thermodynamic Lyapunov function (12.31) generalizing (12.151) to 


FUc,C2X3,---)=F{c{.)) + Y, / ifh(0(x)||C;(x))M,(dx). 


(4.1) 


i>2 


SX3 


The condition of extremality (12.161) extends to 

- [ [ uj{x,y)Xx,dy)Midx)+ [ 

Jx dc{x) Jsx Jsxy 2 




= 0 (4.2) 


and the equilibrium quantities (12.191) become 


Ci(x;c) = c;(x)^ 




(4.3) 


where g is the inverse function to h'. Recall that h was assumed convex on R+ and hence h' is 
an increasing function (0, cxo) —)■ (a, h) with some (finite or infinite) interval (a, h). Hence g is an 
increasing function on (a, b). Let 

G{x) = j g{y)dy 
dyo 


with some yo G [a, b]. Then (7 is a concave function. In particular, for the Burg relative entropy 
(12.41) . h{x) = —\nx and h'{x{) = —1/x is self-inverse, so that g{y) = —1/y. Formula (14.31) 
become 


Ci(x;c) 


c;(x) 

Eiic*(x,)/c(x,)’ 


(4.4) 


and one can choose G{x) = — ln(x). 

Dynamics equations (12.291) or (12.301) remain the same, though of course with ( of form (14.31) 
rather than (12.191) . Consequently the evolution of the free energy remains the same, that is (13.11) 
or (13.21) . The only thing needed a modification is the function 9. Let us restrict the discussion to 
evolution (12.301) and (13.11) only (that is, with restriction (12.241) ) and define 




(4.5) 


so that 


d'W = = f f M{dz)C{z)iy{z-,dy) 

Jsx>2 dsx>2 


6F 

6c{z) 


vMy)/ 


X ^ < A 


6F 

6c{z] 


+ (1-A) 


5F Y 

My)/ 
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Then we get again 0"{X) > 0 and (13.51) . Again condition 6^(0) = 6^(1) turns out to be sufficient 
for the decrease of the free energy by the evolution, and we finally conclude that if the linear 
combinations of the functions 



for all continuous functions c are dense in the space of continuous functions on SX> 2 , condition 
0(0) = 0(1) is equivalent to (13.101) . i.e. to the complex balance condition. 

5. Diffusion approximation 

5.1. Binary mechanisms of diffusion 

Let us consider a lattice hTX^, /i > 0, in R” equipped with the standard basis ei, • • • , e^. To 
each cell or site x = /i(ji, • • • , jn) there is attached a locally compact state space V specifying the 
possible types of particles. Fixing some measure M(dv) in V we can speak about the concentration 
c{x, v) of particles of type v at the site x. The concentration of pairs will be considered with respect 
to the product measure on 

By N{h, x) let us denote the set of neighboring cells to x, that is 


N{x, h) = {y = x± hei, i = 1, - ■ ■ , n}. 


We start here with modeling only the movement of the particles around h7/^, when no change 
of type is possible. We shall assume that only particles in neighboring cells can interact and that the 
interaction is pairwise (which is mostly observed in practice). We shall also assume that our lattice 
is homogenous in the sense that all rate constants, equilibria concentrations, etc, do not depend on 
the site. 

There are three natural mechanisms of transitions between any chosen pair of neighboring cells 
{x,y = X + hci), which we shall also denoted J, II ETl : 

Exchange: (u^, —)■ (u^^, w^), that is, particles of type v, w exchange places; 

Clustering: that is, a particle from one cell attracts a particle from 

another one; 

Repulsion: {v\ w^) (u^, w^^), which is the inverse process to clustering. 

In the spirit of our general approach, we shall assume that any pair of particles, before an inter¬ 
action, should form a compound of two particles. Moreover, the interaction between compounds 
is linear and the number of pairs are in fast equilibrium with the concentration of free particles 
according to the rule (12.201) . that is, the concentration C((3^) w), {Vi w)) of pairs in two neighboring 
cells {x,y = X + hci) equals 



(5.1) 


with some equilibrium C^* (not depending on x, y by the assumed homogeneity) and a free energy 
functional F. 
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For the case of perfect free energy 


F{c{.))= 


In 


c(x, v) 



(5.2) 


this turns to the MAL dependence 



(5.3) 


Furthermore, as we assumed M 2 to have a product form, one can take the transition kernels /r 
from (12.361) . that is dy) = M{dy) and hence (12.141) becomes 



(5.4) 



5.2. Exchange 

Let us start with the reaction of exchange. The linear reaction of the concentrations C{{x,v); {y,w)) 
due to the exchange mechanism between the cells {x, y) is described by the equation 


(((a;,^), {y,w)) = k{v,w)[C,{{y,v), {x,w)) - C,{{x,v), (|/,w))]. 


Here the rates k{v,w) do not depend on the sites by homogeneity, but it can depend on the order of 
the arguments v, w. The r.h.s. of this equation describes the flux of particles along the edge (x, y), 
or, having in mind another equivalent visual picture, through the border of the cells centered at x 
and y. 

Assuming only the exchange mechanism in the system and the MAL condition (15.31) it follows 

that 


c{x,v)='Y] / [C{{x,v),{x + hei,w)) + C{{x,v),{x - hei,w))]M{dw) 
i=i Jy 



X [c(x + /iCj, n)c(x, w) — c(x, v)c(x + hci, w) + c(x — hci, v)c{x, w) — c{x, v)c{x — hci, te)]. 


Introducing the normalized rates 
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and expanding the funetions c in Taylor series up to the seeond order we obtain in the first nontrivial 
approximation 


c{x,v) = / (j){v,w)[Ac{x,v)c{x,w) — Ac{x,w)c{x,v)]M{dw), 


(5.5) 


where the Laplaeian A aets on the first variable of c{x, v). Allowing additionally the evolution of 
free partieles aeeording to the simplest linear dynamies 


C{X,V = 


y / k(v)lc(y,v) — c(x,v)]M(dw) 
dv 


yGN{h,x) ' 

yields in the first approximation the dynamies 

c{x,v) = h‘^k{v)Ac{v, x) + / (p{v,w)[Ac{x,v)c{x,w) — Ac{x,w)c{x,v)]M{dw), (5.6) 


whieh ean be equivalently written in the form 

c{x,v) = h^k{v)diYVc{x,v) + h^ f (j){v,w)dw[Wc{x,v)c{x,w) - Vc{x,w)c{x,v)]M{dw), 

(5.7) 

(with derivations aeting on the first variable of c). 

To get a proper limiting equation one has to assume, of eourse, that k and 0 seale appropriately 
with h, so that the limits 

K{v) = \im h‘^k{v), ^(v,w) = \im h‘^(f>{v,w) 

h^O h^O 

exist, in whieh ease the limiting equation takes the form 

c{x,v) = K{v)div'Vc{x,v) + j ^{v,w)diY\^c{x,v)c{x,w)— Vc{x,w)c{x,v)]M{dw). (5.8) 

The same remark eoneerns all limiting equations below. 

For a more general free energy F(c(.)) the evolution beeomes 


c(a:,u)=^ I k{v,w)(*{v,w)M{dw) 

6F 6F 


i=l 


'V 


X 


exp 


+ 


6c{x + hei,v) 6c{x,w) 


exp 


6F 


+ 


6F 


5c{x — hci^v) 6c{x,w) 


exp 


— exp 


6F 


+ 


6F 


6c{x + hei,w) 6c{x,v) 


6F 


+ 


6F 


6c{x — hei,w) 6c(x,v) j 
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Expanding the variational derivatives in Taylor series up to the second order yields 
5F 5F 


exp 


+ 


5c{x + hci^v) 6c{x,w) 
= exp 


6F 6F ] { a 5F l.d'^dF 

+ — -r (■ exp h—— -^ + -h^ 


Sc{x,v) 6c{x,w) 


= exp 


6F 


+ 


6F 


6c{x,v) dc(x,w) 


1 F h 


dxi Sc(x, v) 2 dx"^ 6c{x, v) 

d 6F 1 2 d‘^ 5F 1 2 ( d 6F 

dxi (5c(x, u) 2 dxj 5c{x, v) 2 6c{x, v) 


and similar with other terms. Thus one sees that zero-order and first order terms again cancel, and 
the second order terms yield the equation 


c{x,v) =h? / /c(u, tu) exp 


5F 


+ 


5F 


6c(x,v) 6c(x,w) 


M{dw) 


X A- 


5F 


6F 


dc{x,v) 6c{x,w) 

which can also be written in the divergence form 



2 


1 

Sc{x, v) 


(5c(x, w) 



(5.9) 


c{x,v) =div J (j){v,w)M{dw 

5F 


X 


exp 


(5c(x, w) 


V exp 


5F 


5c(x, v) 


exp <; j- -r \ V exp 

oc(x, v) 


SF 


6c(x, w) 


(5.10) 


with 


(j){v,w) = h‘^k{v,w)(*{v,w). 


Let us calculate the evolution of the thermodynamic Lyapunov function F(c(.)) along the evo¬ 
lution (15.91) . We shall consider the unbounded lattice hZ^ and its limit (alternatively, one can 
work with finite volume assuming appropriate boundary conditions, say periodic). We have 




6F 


-c{x, v) dxM{dv). 


/Rd Jv 5c{x,v) 

Substituting (15.101) and using the symmetry with respect to the integration variable v, w we get 


1 


2 7^2 [(5c(x,u) 


6F 


(p{v,w) 


6F 


X div 


exp 


SF 


Sc{x, w) 


V exp 


5c(x, w) 
SF 


(p{w,v) 


Sc(x, v) 


dxM{dv)M{dw) 
SF 


exp 


Sc{x, v) 


V exp 


SF 


Sc{x, w) 
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or, integrating by parts in x, 

1 r r 


^*‘*'*' 2 LjyJ [SC{X,V) 


(p[V,W) — — - :-(p[W,V) 


X V 


6F 


6F 


exp 


Sc(x,v) 6c{x,w) 

Hence, if the detailed balance condition 

(j){v, w) = v) 

holds (or equivalently k{v, w) = k{w, v)), then 


Sc(x, w) 

6F 6F 

+ 


dxM{dv)M{dw) 


dc(x, v) Sc(x, w) j 


^ jRd-y^V^ 


SF 


6F 


Sc(x,v) Sc(x,w) 


dxM(dv)M(dw) 


X exp 


SF SF 

+ 


Sc(x, v) Sc(x, w) j 

which is clearly non-positive. 


(5.11) 


(5.12) 


5.3. Repulsion and attraction 

Let us turn to attraction - repulsion interactions. Introducing the rate constants katr{v, w), describ¬ 
ing the process that pushes a particle n to a neighboring particle w, and krep, describing the process 
with a which a particle w can kick out a particle v (siting at the same site as w) to a neighbor¬ 
ing site, we can write the following linear evolution of the concentrations C((^) ^); illy tu)) due to 
attraction -repulsion mechanism between the cells (x, y): 

Ciixy v), (y, w)) =krepiw, v)Ciix, v), (x, w)) + krepiv, w)Ciiy, v), (y, w)) 

- ikatriv,w) -f katriw,v))Ciix,v), iy,w)), 

C((x, x), (x, w)) =katriw, v)Ci{x, v), (?/, w)) + fcatr(v, w)C{{y, x), (x, w)) 

- (krepiv, W) + krepiw,v))Ciix,v), (x,tx)). 

It is worth noting that krepiv, w) and katriv,w) need not be symmetric functions of v,w. Even 
more so, there are natural situations with, say, katriv,w) > 0 and katriw,v) = 0, which means 
that X is a mobile particle and w is not. 

As now we shall have to take into accounts the compounds of particles sitting on the same site, 
(15.41) generalizes to 


c(x,x) = C(a;,w)+ / [C((a^,^'), ix+hei,w))+C,iix,v), ix-hei,w))]M{dw) 

(5.13) 


'V 


i=l 
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Moreover, fast equilibrium eondition (15.11) should be supplemented by the eondition 

Q[[x,v),[x,w)]c) = Q (u, w) exp ^ + 


5c{x, v) 6c{x, w) j ’ 


(5.14) 


with some C*(^, w) that ean be different from (*(v, w), whieh in the ease of the perfeet free energy 
turns to the MAL dependenee 


A// ^ ^ ^ ^ ^ 

C((x, v), {x, w)-,c)= c{x, v)c{x, w). 

c*{v)c*[w) 


(5.15) 


Thus taking into aeeount only the attraetion-repulsion meehanism, using again for simplieity 
the MAL eondition (15.31) . and introdueing the normalized rates 

<l)atr{v,w) = <l)rep{v,w) = krep{v, w), 

the evolution of the eoneentrations beeomes 


n « 

/ [(j)atr{w,v)c{x,v)c{x + hei,w) + (j)atr{v,w)c{x + hei,v)c{x,w) 
i=l Jv 

— {(l>rep{.W, v) + (j)rep{,V, w))c{x, v)c{x, w)]M{dw) 

n « 

+ y' / [(t)atr{w,v)c{x,v)c{x - hei,w) ^ (f)atr{v,w)c{x - hei,v)c{x,w) 

i=i Jy 

- {(j)rep{w,v) + (j)rep{v,w))c{x,v)c{x,w)]M{dw) 


+ y / [(l)rep{w^v)c{x,v)c{x,w)-\-(l)rep{v,w)c{x-\-hei,v)c{x-\-hei^w) 

i=i Jy 

- (0atr(ty, W) + (t>atr{v, w))c{x, v)c{x + /iCj, w)]M{dw) 

n p 

+ y / [(j)rep{w,v)c{x,v)c{x,w) + (j)rep{v,w)c{x — hei,v)c{x — hei,w) 

i=i Jy 

- {(l>atr{w, v) + (patriv, w))c{x, v)c{x - kCi, w)]M (dw). 

Expanding the funetions c in Taylor series, we see that the terms of zero-order and first-order 
in h eaneel. Expanding up to the seeond order we obtain the equation 


f r (9^c d^c 

--h^ yj J 4 >atr{w, v)c{x, v)-^{x, w) + (patriv, w)c{x, w)^;^{x, v) 


dxj 


M{dw) 


n « 

I 

i=i 


4^rep(d^ 1 




d^c 


dc 


dc 


Q ^2 v) + 2^(x, v)-^{x, w) 


M{dw), 
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or in concise notations 


c{x, v) =h^ J [(j)atr{v, w)c{x, w)Ac{x, v) + (j)atr{u!, v)c{x, v)Ac{x, w)]M{dw) 
+ h'^ (prepiv,w)A[c{x,w)c{x,v)]M(dw). 


(5.16) 


The ’repulsion’ part (with vanishing (patr) of this equation ean also be written in the divergenee 
form: 

c{x,v) = h?dS.'v J (j)rep{v,w)'V[c{x,v)c{x,w)]M{dw). (5.17) 

Generalizing, as above for the exchange meehanism, to more general free energy F(c(.)), 
equation (15.161) generalizes to 


C[X, V) = 


\ ^ ^ , ^A 

(patriv, U!)A— -- + 0atr(w, V)A 


6c(x, v) 


+ / 0rep(^^,^^^)Aexp 


5F 


+ 


5c(x, w 
5F 


exp 


5F SF ) 


+ 


5c(x, v) 6c(x, w) j 


6c{x,v) Sc{x,w) 


M{dw). 


M{dw) 


(5.18) 


where 


0a4r(n, W) = h‘^katr{v, w)C{v, w), (j)rep{v, w) = h‘^krep{v, w)C{v, w). 


Similarly to the ealeulations with exehange meehanism above, we find the following law of the 
evolution of F due to the repulsion meehanism (15.181) (taking vanishing (patr in (15.181) 1: 


[Sc{x,v) 


0rep(n) m) T _ - \Yrep\^y^) 


X V 


6F 


6F 


6c{x, w) 

6F 6F 

exp <; - 7 + 


dxM{dv)M{dw) 


6c(x, v) 6c(x, in) J ( 6c{x, v) Sc(x, w) 

Henee, if the detailed balance condition 

4^rep(^tJ^ w'^ (prepi^tV^ v'^ 

holds (or equivalently krep{v, w) = krep{w, v)), then 

1 f 


(5.19) 


F{C{.))=-- Prep{v,w) 


^ , SF SF 

V ( -r + 


Sc{x,v) Sc{x,w) 


dxM{dw)M{dv) 


SF SF 

X exp <; - r + 


(5.20) 


Sc{x, v) Sc{x, w) j ’ 

whieh is elearly non-positive. 
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5.4. Diffusion combined with other reactions 


Suppose that on the sites of the lattice the particles can react according to (12.301) . though only 
pairs of particles can interact producing only two or three particles. Suppose also the free energy 
is perfect leading to MAL with all equilibrium concentration normalized to unity and that the 
simplest product measure M{dv)M{dw) on can be used to measure the concentration of pairs. 
Then the total dynamics comprising diffusion along the spatial variable (including one-particle 
diffusion, exchange and repulsion-attraction mechanism) and reactions on the sites becomes 


c{x,v) = K{v)di'v'Vc{x,v) + L*c{x,v) 

-|-div / ^{v,w)['\/c{x,v)c{x,w) —'Vc{x,w)c{x,v)]M{dw) 


'V 


div / ^rep{v,w)['\/c{x,v)c{x,w) + Vc{x,w)c{x,v)]M{dw) 


'V 


[^atr{v, w)c{x, w)Ac{x, v) + ^atr{w, v)c{x, v)Ac{x, w)]M{dw) 


'v k=2■ 


sv*^ 


c{x, Uj)u{v, w, dui ■ ■ ■ Uk) — c{x, v)c{x, w)v{v, w, dui ■ 

.J=l 


■duk) 
(5.21) 


where all differentiations act on the x variable and L* acts on the second variable. 

Let us stress that the mathematical difficulties in rigorous study of this type of equations in 
general are enormous. In particular, this type includes the full classical Boltzmann equation, for 
which the well-posedness is a well known open problem. 

As a simple interesting example let us describe the case of only two types of particles, V = 
{A, B}, such that the particles of the second type B are immobile (in particular, there is no ex¬ 
change) and act only as catalysis for the branching of A. If the death rate of A is the cor¬ 
responding evolution of the concentration of A (the concentration of B does not evolve in time) 
becomes 

Ca{x) =KAca{x) - ^dCA{x) + ^repdA[cB{x)AcA{x) + Ca{x)Acb{x)] 

(5 22) 

+ ^atrCB{x)AcA{x) +'^Cb{x)ca{x){c^~^{x) - l)z/fc. 

k=2 

Equations of that type are actively studied now in econophysics as models for economic and 
biological growth, the solutions having quite peculiar properties, see e.g. HTll . 

6. Conclusion 

We studied the Michaelis-Menten-Stueckelberg limit (Figure [3]) and found the general form of a 
nonlinear evolutions describing transformations of particles in this limit which combines QSS and 
QE assumptions about transformations of intermediates. 
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The resulting evolution ean be eonsidered as a far reaehing extension to arbitrary state spaees 
of the theory developed by Miehaelis and Menten for the simple enzyme kinetie and by Stueekel- 
berg for Boltzmann’s gas with eollisions. It is developed both for pure jump underlying proeesses 
and for their diffusive limits. It is shown that the eorresponding (generalized) free energy mono- 
tonieally deereases whenever the evolution satisfies either the detailed balanee eondition or more 
generally a eomplex balanee (or eyelie balanee) eondition. The eomplex balanee eonditions fol¬ 
lows from the Markov mierokineties in the Miehaelis-Menten-Stueekelberg limit. 


7. Appendix 

7.1. On pure-jump Markov processes 

Let A be a loeally eompaet metrie spaee. A generator of an arbitrary pure-jump Markov proeess 
(Markov ehain) on X has the form 

Lf{x) = J{f{y) - f{x))u{x, dy) (7.1) 

with a stoehastie kernel u. The dual operator on measures is 

L*y{dx) = J[^{y, dx)y{dy) — iy(x, dy)y{dx)], (7.2) 

so that the evolution of the distributions of the Markov proeess speeified by L is 

/i = L*fi. 

Let a Radon measure M{dx) (i.e. a Borel measure with all eompaet sets having a finite measure) 
be ehosen on A. we say that a bounded measure /r has the eoneentration or the density-funetion 
c G L^{M) if is absolutely eontinuous with respeet to M with the Radon-Nikodyme derivative 
being c, that is 

/ y{dx) = / c{x)M{dx) 

Jv Jv 

for any Borel set V. In order to be able to restriet the evolution /i = L*y on measures with the 
densities, we have to make the following assumption: 

The projeetion of the measure v{y, dx)M{dy) on x, that is the measure v{y, dx)M{dy) 
on A, is absolutely eontinuous with respeet to M or equivalently (by the disintegration of measure 
theory) there exists a stoehastie kernel i){x, dy) sueh that 

h'{y,dx)M{dy) = i>{x,dy)M{dx). (7.3) 


If this is the ease, 

L*[c{x)M{dx)] 


J [c{y)u{x, dy)M{dx) — u{x, dy)c{x)M{dx)], 


(7.4) 
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and the evolution equation /i = L*n in terms of the eoneentrations beeomes 

c{x)= [c{y)i>{x,dy) - c{x)u{x,dy)]. (7.5) 

Jy&X 

Remark 1. The dual to (17.41) equation on functions is 


fix)=[ [f{y) - f{x)Mx,dy). (7.6) 

Jyex 

Its well-posedness (implying the well-posedness for (17.51) ) is investigated under rather general 
conditions on possibly unbounded v in HW . and / l2^ . 

More generally, if we have n loeally compaet metric spaces Xj, j = 1, • • • , fc, a generator of 
an arbitrary pure-jump Markov process on the disjoint union of these Xj has the form 


k „ 

iLf)j{xj) = / ifiivi) - fj{xj))j^j^i{xj,dyi) 

1=1 d 


(7.7) 


with some stochastic kernels Uj=,i. The dual operator on measures becomes 

k „ 

(L*p)j(dxj) = / [ui=,j(yu dxj)p(dyi) - Uj=,i(xj, dyi)p(dxj)]. (7.8) 

1=1 d 

Extending (17.31) we assume that 


ui^jiyi, dxj)M(dyi) = Uj^fxj, dyi)M(dxj). (7.9) 

In this case the evolution of the distributions /i = L*p with p= {pi, - ■ ■ , Pk) can be restricted 
to the concentrations yielding the evolution 


k „ 

Cjixj) = Y^ / [ci(yi)i)j^i{xj,dyi) - Cj{xj)uj^i(xj,dyi)]. (7.10) 

1=1 


In the simplest case when all Uj^^fxj, dyi) have densities Uj^fxj^yi) with respect to Mi, (17.101) 
turns to 

k „ 

= X] / [ci{yi)i^i^j{yi,Xj) - Cj{xj)uj^i(xj,yi)]M(dyi). (7.11) 

1=1 

Of course evolution (17.101) can be considered as a particular case of (17.51) if X is taken to be the 
disjoint union of spaces Xj. 

Recall now that the concentration c*{x) is called an equilibrium for system (17.51) . if 



c*(x)i'(x, dy)] = 0. 


(7.12) 
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If this is the case, and assuming c*(x) > 0 everywhere, equation (17.51) rewrites equivalently as 

c(y) c(x) 


c{x = 


c*(y) 


/yGX 


c*(y) c*(x) 


u{x,dy). 


(7.13) 


For a convex smooth function h{x) let us introduce the ’generalized entropy’ function 


Hh(,c\\c*) = / c*{x)h 


M{dx) 
c*{x) J ^ ’ 


(7.14) 


Assuming that c evolves according to (17.131) and that all integrals below are well defined, it follows 


that 


d 

dt 


H,(c\\c* 



c\y) 


c{.y) 

c*{y) 


c{x) 

c*{x) 


i>{x, dy)M{dx). 


(7.15) 


Generalizing the concepts from the theory of Markov chains let us introduce the graph (X, E) 
associated with evolution (17.41) such that the set of vertices X coincides with the state space X and 
the edge {x —)■ y) exists if the point y belongs to the support of the measure ^{x, dy). As usual, 
we say that the finite sequence (|/o, 2/i, • • • , 2 /^) is a path in this graph joining y^ and yk if the edges 
{yj-i —)■ yj) exist for all j = 1, • • • , fc; and that the graph is strongly connected if for any pair of 
points {yo, y) there exist paths joining y^ and y. 

The following result is the extension of the Morimoto H-theorem of finite state-space Markov 
chains to the continuous state-space: 


Proposition 2. Under evolution (17.131) . and assuming c*{x) > 0 everywhere, 


dHh{c\\c*) 

dt 


< 0 . 


(7.16) 


Moreover, if the measure M(dx) has the full support and the graph (X, E) introduced above is 
strongly connected, then the equality in (17.161) holds if and only if the ratio c(x)/c*(x) is a constant. 


Proof As it follows from (17.21) . / L*p{dx) = 0 for all p. In terms of equation (17.131) this rewrites 
as 


0 


c{x)M{dx) 


\y)u{x,dy) 


c{y) 

c*{y) 


c{x) 

c*{x) 


M{dx) 


= / c*{y)u{y,dx) 


c{y) c(x) 


c*(y) c*(x) 


M(dy) 


(7.17) 


for any c(x). Consequently, for any function / (such that the integral below is well defined). 


[ c*{y)u{y, dx) [f{y) - f{x)] M{dy) = 0. 
Jx 

This identity allows one to rewrite (17.151) as 

^^Hh{c\\c*) = J J c*{y)v{y,dx)M{dy) 


(7.18) 
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X 


h 


c{x) 


h 


c{y) 


h' 


/ f c(a:) \ f c{y) c(x) 


c*(x) J \c*{y) c*{x) J \ 


(7.19) 


c*(a;) J \c*{y) 

implying (17.161) by the convexity of h. 

Finally, assuming M has full support, it follow that the equality in (17.161) holds if and only if 


u{y,dx) 


h 


c{x) 

c*{x) 




c*{y) 


c*{x) J \c*{y) c*{x) 


= 0 


for all y. Hence by convexity, for all x from the support of v{y, .). The final conclusion 


c*{x) c*{y) 

follows from the assumed connectivity of (X, E). 


□ 


7.2. Linking the concentration of particles and of compounds 


Proposition 3. The kernels fi in (12.1 II) can be chosen in such a way that ifCki^i, • • • , Xk) = Ck{^) 
is the concentration of the compounds x of size k, the concentration of particles involved in these 
compounds equals (12.131) . that is 


c{x)= / Ck{x,X2,-■ ■ ,Xk)pk{x,dX2-■ -dXk). 


(7.20) 


Proof. Let firstly k = 2. The arbitrary measure M 2 on SX'^ can be given by the pair of measures 
and (the subscripts d and nd stand for diagonal and non-diagonal parts), where is a 
measure on the diagonal D = {{x,x) : x G X} and is a symmetric measure on X^ \ D, so 
that, for a symmetric function /, 


/ f{x,y)M 2 {dxdy) = - f{x,y)M^‘^{dxdy)+ f{x,x)M^{dx). 

'SX2 ^ 7x2 7 


(7.21) 


Assuming that M 2 has absolutely continuous (with respect to M) projections on X means that 
there exist a kernel p'^^{x, dy) with /i”'^(x, {x}) = 0 and a function a;(x) such that 

M"‘^{dxdy) = M{dx)p"‘^{x,dy), M^{dx) = uj{x)M{dx). 

Then clearly (17.211) becomes 

1 


/ f{x,y)M2{dxdy) = 

'sx^ Jx ^ 


f{x,y)p2{x,dy) 


M\dx) 


(7.22) 


with 


/i 2 (a:, dy) = p^'^{x, dy) + 2a;(x)5(x - y). 

Moreover, the amount of particles in a neighborhood dx of a point x entering the compounds is 

C{x,y)M^'^{dxdy) + 2 f ({x,x)M‘^{dx) 

f dx J X J dx 



= / M{dx) 

J dx 


C,{x,y)tE (x,7|/) + 2C(x,x)a;(x) 


' X 
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(a particle at x is used twiee in the compound ({x, x), hence the eoeffieient 2 at the seeond term). 
Hence the concentration, whieh is the density with respect to M{dx) is 


c{x) = / C{x,y)iJ.'^'^{x,dy) + 2C{x,x)uj{x) 
Jx 


C{x,y)y2{x,dy), 


as required. 

Now let k = 3. Then an arbitrary measure M 3 on SX^ ean be given by the triple M*^, M”'^ 
and M*”^ where M'^ is a measure on the diagonal = {{x,x,x) : x E X}, M"'^ is a symmetric 
measure on X^ \ where 


= {{xi,X2,X3) : 3i,j : Xi = Xj], 


and M*'^* is a measure on X'^ (not neeessarily symmetrie, that eounts the triples (x, x, y) with 
y ^ x) so that for a symmetrie funetion /, 



f{xi,X2, x^)M^{dxidx2dx^) 


1 

6 



/(xi, X2-, x^)M'^'^{dxidx2dxz) 


+ / f{x,x,x)M‘^{dx) + 

Jx 


'{XxX)\D 


f{x,x,y)M^^ {dxdy) 


(7.23) 


Assuming that M 2 has absolutely eontinuous (with respeet to M) projeetions on X implies that all 
three measures above have this property and the proof of the statement ean be performed separately 
for each of them. For M"^'^ and M'^ it is literally the same as for the case k = 2. Let us eonsider 
a more subtle ease of the measure Denoting by yu and ^21 the kernels arising from the 

projeetions of M*”* on the first and the seeond eoordinate (note that they are not symmetrie, as the 
first eoordinate deseribes the pairs of identieal partieles), we have 

M’‘"‘^{dxdy) = M{dx)y 21 (x, dy) = M(dy)y^iy, dx) 


and therefore also 

M'-'^^dxdy) = ‘^M{dx)fi 2 i{,x,dy) + ^M{dy)yi 2 {y,dx). 

Consequently, defining the kernel 

y{x,dydz) = 2y2i{x, dy)6{z - x) + yi 2 {x, dy)6{z - y), (7.24) 


allows one to write 


IXxX\D 


f{x,x,y)M^^ (dxdy) = 


'X 


'SX 2 


f(,x,y,z)y{x,dydz) 


M{dx) 


Moreover, the amount of particles in a neighborhood dx of a point x entering the compounds 
containing precisely two identieal partieles equals 



C{x,y,y)M^^ {dydx) + 2 


dx JX 



C{x,x,y)M^^\dxdy) 


dx JX 
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= / M{dx) 


' dx 


C,{x,y,y)iii2{x,dy) + 2 / C,{x,x,y)y2i{x,dy) 


IX 


Hence the concentration, which is the density with respect to M{dx) is 


c{x) 



C{x,y,z)y{x,dy dz), 


as required. 

Larger k are analyzed similarly, but requires understanding of the structure of measures on 
SX^ discussed below. □ 


Recall that a partition of a natural number k is defined as its representation as a sum of non¬ 
vanishing terms (with the order of terms irrelevant), i.e. as 


k — Ni -f 2 N 2 + ■ ■ ■ + jXj 


(7.25) 


with a j > 0, where Ni is the number of terms in the sum that equal 1. Graphically these partitions 
are described by the so-called Young schemes. For a partition (or a Young scheme) (17.251) let us 
defined the extended diagonal as a subset of the product such that at least 

two of the coordinates (xi, X 2 , • • • , xn^^+.-.+Nj) coincide. The following fact is then more or less 
straightforward. 

Proposition 4. An arbitrary Borel measure Ms on SX^ can be uniquely specified by a collection 
of measures on X^^^ \ which are symmetric for permutations inside the 

group of arguments in each X^'- and which are parametrized by all partitions (17.251) . so that for a 
symmetric function f on X^ 


[ f{x)Ms{dx)= Y. 

JSXl^ .r ___ JVl! • • • 




■Nj\ 


M^^’-X(^dxi ■ ■■dxNi+-+Ni) 


^ f i,X\ 5***5 Xj\[-1 5***5 X -* * * 5 -* * * ) 


Itimes 


(7.26) 


(the arguments coincide in each group entering the partition), the sum being over all partitions 
(17.251) of k. If additionally, the projection of Ms on X is absolutely continuous with respect to 
a measure M(dx), that is each measure is absolutely continuous with respect to each 

arguments, then it can be presented in Ni + ■ ■ ■ + Nj equivalent forms: 

M^^’-X(dxi ■ ■ ■dXNi+.-.+Nj) = M{dXNi+-+Ni_i+m) 

Pi (X]\f-^-\ -dxi • • • dXjViH- l-Ni_i+m ' ' ' dx -(7.27) 

where Xp denotes, as usual, the absence ofxp in the sequence of arguments, pi are some stochastic 
kernels and m G {1, • • • iV)}, or more symmetrically as 

j ^ Ni 

M^^’-Xj(dxi ■ ■ ■ dXNi+-+Nj) = X] ^ X] M{dXNi+-+Ni_i+m) 

1=1 m=l 
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^ {^Ni+---+Ni_i+m,dXi- ■ ■ dX]Sf^+...+Ni_i+m- ■ ■ dXNi+---+Nj)- (7.28) 

The numerators I in (17.281) refleet the number of identieal partieles entering a eompound, thus 
presenting the analogs of stoiehiometrie eoeffieients. 
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